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Abstract 

In this paper we introduce a new method to study the influence of thermal 
conduction and viscous processes in anisotropic gravitational collapse. To this 
end we employ the HJR method to solve the Einstein equations. The Maxwell- 
Cattaneo type transport equations are used to find the temperature and bulk 
and shear viscous pressures. Under some conditions Maxwell-Cattaneo trans- 
port equations comply with relativistic causality. Thus, it is advisable to 
use them instead of Eckart transport equations. In the inner layers of the 
star the temperature ceases to be sensitive to the boundary condition. This 
behavior, which can be explained in terms of the Eddington approximation, 
allows us to find the thickness of the neutrinosphere. The dynamics of col- 
lapsing dense stars is deeply influenced by the neutrino emission/absorption 
processes. These cool the star and drive it to a new equilibrium state. There- 
fore, the calculation of transport coefficients is based on these processes. 
04.25.Dm, 05.70.Ln, 97.60.-s 
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I. INTRODUCTION 



The physical study of relativistic stars rests usually on the assumption of local isotropy. 
However, at least for high densities, theoretical evidence suggests that this may not be a very 
accurate approximation 0,0]. Then, it seems fitting to adopt anisotropic models to carry 
out a more detailed analysis of the internal structure of neutron stars, their upper-mass limit 
and their collapse dynamics. Because of the anisotropy, shear viscous stresses arises. The 
interchange of momentum between the different layers of the star causes viscous processes. 
It unavoidably comes from the interaction between particles and from the matter-radiation 
interaction as well. Therefore, the presence of radiation flows during the collapse reveals the 
existence of anisotropy in the system. When radiation and matter interact strongly (i.e., 
diffusion) shear and bulk viscous pressure arise. On the opposite limit (free streaming) one 
has anisotropy associated to the direction of propagation of the radiation. Because of this, 
it is suitable to introduce an anisotropic state equation if one wishes to analyze systems in 
which the flow of radiation non negligibly influences the dynamics of the collapse It is 

then important to study the anisotropy associated to the radiation flow during the collapse. 

We shall adopt the semi numerical method of Herrera et al. [^] - HJR henceforth - to solve 
the Einstein equations. Furthemore, to find the temperature and bulk and shear viscous 
pressure it will be necessary to introduce a set of transport equations. To solve them we will 
adopt the explicit expressions for the transport coefficients of heat conductivity and bulk 
and shear viscosity of an interacting mixture of matter and radiation found by Weinberg . 
This requires to determine which reactions are more relevant in the momentum transport 
between the different star layers. The large neutrino mean free path implies high transport 
coefficients. Thus, viscous processes induced by interactions between neutrinos and matter 
are much more important than those induced by other scattering processes. The cooling by 
absorption and emission of neutrinos drive the star to the equilibrium. Among the sources 
of neutrino opacities it seems that the scattering with electrons (e" + z/ — + z/) and 
nucleon absorption + n + p) may account for this process pl-pl. 
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The temperature is a key quantity to decide which processes can take place during the 
collapse. Moreover it is indispensable to calculate the bulk and shear viscous pressure. 
Unluckily the solution of the Einstein equations does not provide any information about it. 
This is why a set of transport equations must be adopted. The standard Eckart's theory of 
irreversible processes []T2|-|I1] exhibits undesirable effects (i.e., it predicts an infinite speed 



of propagation for thermal and viscous signals and unstable equilibrium states). It is then 
necessary to resort to another thermodynamic theory of irreversible processes that does 
not present this anomalous behavior. The extended irreversible thermodynamics theory 



(EIT) |ir5| , p!6| seems to be a good candidate to replace the old Eckart theory and it has 



been used in cosmological problems with good results [|T7| -[T9[|. Essentially the EIT rests on 
two hypotheses: (1) The dissipative flows (heat flow and viscous pressures) are considered 
as independent variables, hence, the entropy function depends not only on the classical 
variables (particle number and energy density) but on these dissipative flows as well. (2) 
At equilibrium state, the entropy function is a maximum. Moreover, its flow depends on all 
dissipative flows and its production rate is semipositive definite. As a consequence the heat 
flow Q^, the bulk viscous pressure 11 and the traceless viscous tensor ti^^ obey the evolution 
equations [^] 



TQ Q" K + Q^ = xh^" 



rnt/"n„ + n = -ct/;^ + aoCQf;. 

and 

respectively, where rg, rn and r^r are the relaxation time of thermal, bulk viscous and shear 
viscous signals respectively. T stands for the temperature, x fo^' the thermal conductivity 
coefficient, and Q and rj for the bulk and shear viscous coefficients respectively. The param- 
eters ao and ai are connected with the characteristics of the fluid under study. The spatial 
projector tensor is designed by h^y = g^y - U^Uy, and a^^ = Ui^^-y)- U{^i U^) - hf,yU^^/3 
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is the shear tensor. Here, the upper dot denotes Aai3...= U'''Aaf3...-\- It can be seen that 
if some reasonable assumptions are met, these equations preserves the causahty condition 
and predict stable equilibrium configurations []2T|. However, the transport equations for 
the dissipative fiows predicted by EIT are extremely involved. Nevertheless, in the linear 
approximation the evolution equations for dissipative fiows reduces to 



TQ Q" K + - Xh^" Tu -TU, 
rnt/"H,, + H ^ -(U^^ 

and 



which are the covariant formulation of Maxwell-Cattaneo transport equations ||2^j23[] . One 
may hope that the propagation speed of the dissipative signals is comparable to the sound 
velocity. For this to be true the relaxation time in the transport equations must be nearly 
the radiation mean free time. After calculating the temperature we shall obtain both, the 
bulk and shear viscous pressures, which cannot be directly found from Einstein equations. 
In this way we will establish the infiuence of shear viscous pressure on the anisotropy of the 
system. 

The paper is organized as follows. In section II the Einstein equations, for an anisotropic 
sphere composed by a mixture of radiation and matter, are derived, and we resort to the HJR 
method to solve them. In section HI is we write the Maxwell-Cattaneo transport equations 
within the HJR formalism and find the mean free path of the neutrinos as a function of the 
temperature. This is applied to a specific case in section IV. As initial configuration, we 



shall adopt the anisotropic static solution of Gokhroo and Mehra . This one corresponds 
to an anisotropic fiuid with variable energy density. Lastly, section V summarizes the results 
of this work. 

We adopt metrics of signature —2. The quantities subscripted with 'a' denote that they 
are evaluated at the surface of the star. The subscripts and 1 indicate partial differentiation 



with respect to time (m) and radial coordinate (r), respectively, finally an upper dot on 
scalar quantities means d/du. A hat over a quantity means that this one is measured by a 
Minkowskian observer comoving with the fluid. 



II. FIRST SET OF EQUATIONS: EINSTEIN FIELD EQUATIONS 

A. Interior and exterior metrics 

Our aim is to describe a non-static spherically symmetric fluid distribution. To this end 
we adopt the radiation coordinates [25,^|, then the interior metric takes the form 



ds 



2 _ ^2/3 



V 2 

— du + 2dudr 
r 



- (d9^ + sin^dd<f) , (2.1) 



where u = x° is a timelike coordinate, r = is the null coordinate and 9 = and 
(f) = are the usual angle coordinates. The w-coordinate is the retarded time in flat 
space-time and, therefore, the u-constant surfaces are null cones open to the future. The 
metric functions (3 and V in equation (|2.1| ), are functions of u and r. A function m{u,r) 
can be defined by 

V = e'^'^ir -2m{u,r)), (2.2) 

which is the generalization, inside the distribution, of the "mass aspect" defined by Bondi 
et al In the static limit it coincides with the Schwarzschild mass. 



In order to give a clear physical meaning to the above formulae, we introduce local 
Minkowskian coordinates (t, x, y, z) related to Bondi's radiation coordinates by 

dt = e^sj^du + \l^dr), (2.3) 



dx = e^\j^dr, (2.4) 



dy = rde, (2.5) 
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dz = rsiw 



(2.6) 



and to the Schwarzschild coordinates (T, R, 9, $) by 



T = u+ -dr, 
Jo V 



(2.7) 



R = r Q = e, $ 



(2.8) 



Outside matter the metric is the Vaidya one [^], a particular case of the Bondi metric 
in which j3 = and V = r — 2m. 



B. Stress-energy tensor 

As we mentioned above we consider an anisotropic fluid sphere composed by a material 
medium plus radiation. The material medium travels in the radial direction with a velocity 
uj ma, Minkowski coordinate system. The radiation is treated in the diffusive regime because 
it interacts strongly with matter. Due to this interaction, viscous processes appear. Thus, 



a viscous term must be considered in the stress-energy tensor. We denote this term by 
and, for a local Minkowskian observer comoving with the fluid, it can be written as 



V 



(2.9) 



where 7r^,y denotes the traceless viscous pressure tensor, h^j^y = rj^i, — U^Uu the spatial 
projection tensor and 11 is the bulk viscous pressure. For this comoving observer f/^ complies 
with = 5* , and applying the condition tt^^U^ = together with traceless character of 



TT/iiy, we obtain 



vr 



\ 





TT 










(7r/2) 








[n/2) 



(2.10) 
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Not all anisotropies can be explained in terms of the interaction between matter and ra- 
diation. Therefore, it is necessary to bear in mind that the material part of the stress-energy 
tensor must reflect the anisotropic character of the fluid. For the Minkowskian observer the 



anisotropic material part of the stress-energy tensor T^^ is given by the expression 



T^f, = ipM + Pt)U^U, - PtVf^u + (P - Pt) X^X. 



(2.11) 



where x^l = Pm denotes the material energy density, P the material pressure and Pt = 
P + p the material part of the tangential pressure. Thus, the pressure p refers to the 
anisotropies that cannot be explained as a result of the matter-radiation interaction. 

The last term in the stress-energy tensor accounts for the presence of radiation. In the 
Lagrangean frame (the proper frame) and in the diffusive regime it reads 

/ \ 



T, 



R 



V 



PR 


-Q 








-Q 


V 














V 














V 



(2.12) 



where pr denotes the radiation energy density, Q the heat flow and V the radiation pressure. 

Thus, for a local observer comoving with matter, with radial velocity cj, the stress-energy 
tensor in local Minkowskian coordinates is 



rp rpM I rpR . rpV 



and in virtue of (|2.9|), ( ^.101 ), ( ^.llj ) and ( |2.12| ) it can be written as 



f^, = {p + Pa_)U^U, - P±r],.u + {Pr - P±)Xf.Xu + 2Qi^A) 



(2.13) 



with 



f/^ = (1,0,0,0), 



(2.14) 



= (0,1,0,0), 
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(2.15) 



Q^ = Q(0, -1,0,0), 



(2.16) 



P = PM+PR, (2.17) 
P^ = P + p + n + 7r, (2.18) 



and 



P± = Pr~^TT + p. (2.19) 

The physical variables (energy density p, radial pressure Pr, tangential pressure Pi and heat 
flow Q) are obtained as measured by the mentioned observer, and the effects of gravitation 
are clearly provided through the appropriate transformation to a curvilinear coordinate 
system. 

To study the dynamics of the system (i.e., uj) it is necessary to obtain the stress-energy 
tensor as seen by an observer at rest with respect to the Minkowskian coordinates. To 
this end it is necessary to apply a Lorentz boost in the radial direction to T^^". A further 
transformation allows us to express the stress-energy tensor in curvilinear coordinates (see 
H for details). Thus, applying a local coordinate transformation (|2.3[-|2.6|) and a Lorentz 



boost in the radial direction to ( ^.13| ) we obtain the stress-energy tensor 



T^u ={p + P±)U^U, - P±g^, + (Pr - P±)x^.Xi^ + 2Q(,.t/.), (2.20) 

as measured by an observer using Bondi coordinates with a radial velocity, with respect to 
the matter configuration, —uj. After performing the transformation, expressions ( |2.14| - prT6D 
reads 

and 



2.22 



(2.23) 



Note that Q^'U^ = and Q =^-Q^'Q^. 

Applying these transformations to the traceless viscous tensor (|2.10|) , we get 



TT 



,2/3 



V 



















-TT 













2 / 

Outside matter the stress-energy tensor corresponds to that of a null fluid, i.e. 



where 



-2/3, 



(2.24) 



(2.25) 



(2.26) 



C. Einstein field equations 

Inside matter the Einstein field equations, G^j^^ = SttT^^, can be written as: 

1 ....... X 1 



47rr(r — 2m) 



(^—rrioe + (1 — 2m/r)mi 
1 



mi 

r — 2m 1 — UJ 



{p-2Q + Pr), 



{p + 2ujQ + PrU'), 



l+LU 



- ^ + ^(1 - 2^)(2/3n + - ^) + '^^^^"f^^"^" 
An Sit r r Svrr 

while outside matter, the only non vanishing Einstein equation yields 

2m{u) \ 



Pi 



mo 



-A-nr e 1 



(2.27) 
(2.28) 
(2.29) 
(2.30) 

(2.31) 



To algebraically solve the physical variables present in the above set of field equations 
( [^.27] - PT^ ) we must find out both (3 and m. The HJR method will help us in this task, but 
before it we must impose the junction conditions between exterior and interior metrics. 



D. Junction conditions and surface equations 



The mass function can be expressed as 

rh = r A-nr'^pdr. (2.32) 







It involves an effective energy density given by the right-hand side of ( p. 28 ) 



P = -^{p - Q(l - ^) - Proo), (2.33) 



which in the static hmit reduces to the energy density of the system. 



From ( p.29|) one has 



/3 = /\ ^^^(P - 2Q + Pr)dr, (2.34) 
Ja{u) r — 2m 1 + UJ 



and we may re-write the non-static case as 



rr 27rr'^ 

(3= (~p + p)dr, (2.35) 

Ja{u) r — 2m 



with 



P = -^{-up-Q{l-u;) + Pr) (2.36) 

1 + UJ 

being the effective pressure, which also reduces to the radial pressure in the static limit. 

Matching the Vaidya metric to the Bondi metric at the surface ( r = a) of the fluid dis- 
tribution implies f3a = (3{u,r = a) = with the continuity of the mass function m{u,r) (i.e. 
the continuity of the first fundamental form) and the continuity of the second fundamental 
form leads to 

.rha\ Pa 



1-2— ] ^ " . (2.37) 

J Pa+Pa 



see 



pO| for details). 

From the coordinate transformation (|2.3| ) the velocity of matter in Bondi coordinates 
can be written as 
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dr V uj 
du r 1 — ' 

evaluating the last expression at the surface and comparing it with ( p.37|) it follows 



(2.38) 



Pa — —^aPai 



(2.39) 



or equivalently using ( |2.33| ) and ( |2.36|) we get 



O = P 

=<^a ra) 



(2.40) 



which is a well known result for radiative spheres |31[ . 

To derive the surface equations we introduce five dimensionless functions 

1 



A = 



m(0) 



M = 



m 



u 



u = 



m(0) m(0) 



2M ^ 
F = l Q = , 



f2.411 



where m(0) is the initial mass of the system. Using the functions just defined into (|2.37|) we 
get the first surface equation 



A = F{Q-1). 



(2.42) 



This first surface equation gives the evolution of the radius of the star. 

The second surface equation emerges from the luminosity evaluated at the surface of the 
system. The luminosity as seen by a comoving observer is defined as 



(2.43) 



Evaluating ( ^.27| ) and ( |2.28D at the surface and using the expansion 



mo^ ~ m — ami^, 



the luminosity perceived by an observer at rest at infinity reads 



(2.44) 



L = -M = E{2n - 1)F = AnA^Qa {2Vt - 1) F. 



The function F is related to the boundary redshift Za by 



(2.45) 
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lj^Za = — = F-^'\ (2.46) 

^Tec 

Thus the luminosity as measured by a non comoving observer located on the surface is 

E = L{l + Zaf = -^ = E{2^l-l), (2.47) 
r 



where the term — 1) accounts for the boundary doppler-shift. Using relationship ( p. 45 ) 
together with the first surface equation we obtain the second one as 

which express the evolution of the redshift at the surface. 

The third surface equation is model dependent. For anisotropic fluids the relationship 
(T^fj)a = can be written as 

^ ^ + P ^ +/?j_„-f-(P,,-P)') =0, (2.49) 



\1 — 2m /r y " \r 

where 

It is possible to associate a physical meaning to R±a'- the first term is tied to the hy- 
drodynamic force, the second one to the gravitational force and the third one reflects the 
anisotropic character of the fluid. The latter is negative when the radial pressure is larger 
than the tangential pressure. In this case the fluid distribution may become more compact 
than in the isotropic one. 

Using the expansions (|2.44|) and 



(P + ^)0. ^ [pail - ^a)]o - a{P + p)i,, 



where ( |2.39| ) has been used, we get 



ra 



[I -a) 



F II Pa Pa A p, 

AT^Apa—: -— + Fl^^ + - p. 

il 2A Pa Apa 



r )a 



(2.51) 
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Expression ( p.51|) generalizes the Tolman-Oppenheimer-Volkov equation to the non-static 
radiative anisotropic case. 

Now we are in position to brief summarize the HJR method [^] apphed to the anisotropic 
case. 



E. The HJR method 

This method has been extensively used since last decade to obtain models of non static 
radiating fluids from known static solutions of the Einstein equations - [^,^^ and references 
therein. Its use reduces the problem of the gravitational collapse to the solution of the surface 
equations ( p.42[p3l| ). This procedure simplifies the resolution of the problem because these 
ones are a set of ordinary differential equations that can be solved by Runge-Kutta method. 
The algorithm can be summarized as follows. 

1. Take a static but otherwise arbitrary interior solution of the Einstein equations for a 
spherically symmetric fluid distribution 

Pst = P{r), pst = p{r). 

2. The effective quantities p = p{u,r) and P = P{u,r) must coincide with pst and Pgt 
respectively in the static limit. We assume that the r-dependence in effective quantities 
is the same that in its corresponding static ones. Nevertheless, in non static case the 
junction conditions are different accomplishing expression (|2.39| ). This condition allows 
us to find out the relation between the w-dependence of p = p{u, r) and P = P{u, r). 

3. Introduce p{u, r) and P{u, r), into ( p.32| ) and (|2.35|) to determine m and (3 up to three 
unknown functions of time. 

4. The three surface equations form a system of first order ordinary differential equations, 
by solving it we find the evolution of the radius, A{u), and two unknown functions of 
time. These ones can be related with the w-dependence of p = p{u, r) and P = P{u, r). 
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5. There are four unknown functions of time {A, F, Q and L). Thus, it is necessary to 
impose the evolution of one of them to solve the system of three surface equations. 
Usually the luminosity is taken as an input function because it can be found from 
observational quantities. 

6. Once these three functions are known, it is easy to find rh and f3. Therefore, the interior 
metric is completely defined. 



7. Now, the left-hand side of the Einstein equations (|2.27| - |2730| ) is known. However, 



the right-hand side of these equations contain five unknown quantities (tu, p, Pr, P± 
and Q). Thus, it is necessary to supply another equation to close the system of field 
equations. In the anisotropic static case it can be found a general equation that relates 



the tangential pressure to the mass function, energy density and radial pressure 33 



P.-..4p..(^^)(^). (2.52, 

This expression is usually generalized, in the context of HJR method, to non-static 
cases by substituting for the effective variables the physical quantities PJ5^ 



Now, the Einstein equations, supplemented with (|2.53|) , make up a close system of 
equations and quantities u, p, R,., P± and Q can be found. 

The HJR method starts from a static interior solution of the Einstein equations ( Pgt = 
P{r) and pst = pi'f'))- The effective variables ( p.33| ) and (|2.36| ) must reduce in the static case 



to Pst and Pst respectively. Nevertheless, these quantities have a dependence on the retarded 
time u. We assume that the effective variables have the same r-dependence as the physical 
variables of the static situation have. This can be justified in terms of the characteristic 
times for different processes involved in the collapse scenario. If the hydrostatic time scale 
Trydri ~ ^1 \JGp^ is much shorter than the Kelvin- Helmholtz time scale, Tkh^ then in first 
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approximation the inertial terms in the equation of motion T/f = can be ignored 

L 'f^J r=a 

p5| , page 11]. The Kelvin- Helmholtz phase of the birth of a neutron star last some tens 
of seconds |^6[, whereas for a neutron star of one solar mass and a ten- kilometer radius, 
we obtain Thydr ~ lO^^s. Thus, in this first approximation the r-dependence of physical 
quantities P and p are the same as in the static solution. The assumption that not the 
physical quantities but the effective variables ( 2.33| ) and ( |2.3(i| ) have the same r-dependence 



as the physical variables of the static situation, is a correction to that first approximation. 
Therefore it is expected that the introduction of two functions of time in p and P preserving 
the same r-dependence as in p^t and Pst will yield good results. These two functions of time 
can be related by means of the junction condition ( |2.39|) . 



If the evolution of effective variables where known the functions m and (3 could be found 
and the Einstein equations supplemented by ( p.53| ) would constitute a closed system of 
differential equations. This method will be clarified by means of an example in section IV. 

III. SECOND SET OF EQUATIONS: TRANSPORT EQUATIONS 

We shall use the method described above to find the heat flow, the energy density, the 
radial and tangential pressures and the velocity of the collapse. However, if one wished 
to explicitly find the temperature and viscous pressure, one has to resort to the transport 
equations laid down by EIT. 

Usually classical theory |ll2|-[T^ has been employed as a first approximation to the study 



of gravitational collapse. Nevertheless, this one presents two important disabilities |^ll: (1) 
It predicts an infinite speed for thermal and viscous signals. (2) The equilibrium states turn 
out to be unstable. Therefore, as a further step in the thermodynamic study of gravitational 
collapse, it seems suitable to resort to a theory free of such drawbacks. As we pointed 
out in section I the relativistic EIT may do the job. The Maxwell-Cattaneo transport 



equations are a particular case of those given by EIT -see for instance p2| , p3| . The former 



are decoupled expressions of the EIT transport equations and they may be understood as 
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a first approximation to them. Thus, we consider that the dissipative flows are decoupled 
and the transport equations may approximated by their Maxwell-Cattaneo form. The heat 
conduction equation can be written as 



TQ Q' Ul + ^ xh^" 



(3.1) 



where h^'^ = g'^'^ — U^lf^ is the spatial projection tensor, x is the thermal conductivity coef- 
ficient, T the temperature and tq is the relaxation time for thermal signals. The evolution 
of bulk viscous pressure 11 is given by the expression 

rnf/"n „ + n ~ -Ct/^. (3.2) 

The bulk viscosity coefficient and relaxation time for bulk viscous signals are denoted by ( 
and Tjj respectively. On the other hand the evolution of shear viscous pressure tt^^^ ( |2.24| ) 
can be written in Maxwell-Cattaneo form as, 

Ttt t^^lu +T^^iu - 2770-^1., (3.3) 

where shear tensor h^^Uj^y^a ~ U^Ji^yj?) is denoted by a^^, corresponds to relaxation time 
for shear viscous signals and 77 indicates the shear viscosity coefficient. 
Eckart transport equations are given by 



(3.4) 



n ^ -C^^^ (3.5) 

and 

Ti^u ~ 2ria^u. (3.6) 

Note that the main difference between both formulations is the presence of the relaxation 
time in ( |3.1[ - |3.3| ). As we will show in section IV this discordance is important close to surface. 

To solve this system of partial differential equations it is necessary to adopt an expression 
for transport coefficients x^ V C- Fo^' ^ mixture of matter and radiation these ones are 
given by (see e.g. 0) 
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X 



(3.7) 
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-6TV 



(3.8) 



and 



C = 15r/ 



9p, 



(3.9) 



where r denotes the mean free time of radiation (in our case neutrinos), while the constant 
b takes (for neutrinos) the value 7a/8, where a is radiation constant. 

Expressions ( |3.1[ - |3.3| ) implies a finite propagation speed for viscous and thermal signals. 
Nevertheless, to guarantee relativistic causality it is necessary to restrict the possible values 
of the different relaxation times. If we assume that the propagation speed of thermal and 
viscous signals are of the order of sound speed, then 



T. 



(3.10) 



This means that the relaxation times are somewhat larger than the mean free time of 
radiation. Thus, if r is known, it is possible to use ( |3.1| ) to find the temperature and, 
by (p.2|) and (p.3[) to calculate the viscous pressures. 



A. Determination of r for neutrinos 

Thermal neutrino processes are important during the late stages of collapse of massive 
stars P]. They help to carry thermal energy from the core to outer regions. Thus, we assume 
that neutrinos are principally thermally generated with energies close to ksT. 

A detailed study of the contribution of different processes to the neutrino transport has 
been given by Bruenn [^]. Our aim is to establish an approximate expression to the mean 
free time for neutrinos. Despite the procedure used in this work yields good results, it is 
highly idealized and only approximate. A more accurate expression for r would be desirable, 
but we consider that the approach below describes the process accurately enough. 
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At the typical densities of a neutron star, matter can capture neutrinos. The neutrino 
trapping helps to drive the neutrinos to local equilibrium. The matter opacity to the neu- 
trinos takes place principally by means of electron-neutrino scattering (e~ + u —>■ e~ + u) 
and nucleon absorption {u + n —>■ p + e~). 

As a first approximation the cross-section for nucleon absorption may be taken as 
chapter 18] 

an ~ lO'^^el cm^ (3.11) 

where e^, is the energy of neutrinos in MeV. Then, the neutrino mean free path is given by 

1 10^0 

An = ~ — 5" cm, (3.12) 



where n„ 6 x lO^^p cm~'^ P,|37|, is the neutron number density. 

In high energetic collisions the cross-section for the electron-neutrino scattering can be 
approximated by 

cTe ^ 10-^^e^ cm^. (3.13) 

The electron number density Ue, can be related to by means of the electron fraction 
{jie — Unye)- Heucc, the mean free path is 

1 10^0 

Ae = ^ ^ cm. (3.14) 

ne(Te pe^ye 

The effective mean free time which takes into account the neutrino zig-zag path may be 
written as 

cm. (3.15) 



Assuming that the neutrinos are generated by thermal emission {e^, ~ ksT) we find the 
r dependence on temperature 



r oc 



T-3/2. (3.16) 
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A rigorous treatment of mean free time would comprise the evolution of 3^e- Nevertheless, 
this effort doesn't seem advisable here due to the degree of the approximation adopted and 
that the range of possible values for ye is severely restricted (0.2 < 3^e < 0.3) 0. 

It is convenient to express r in dimensionless form to simplify the numerical treatment 
of equations ( |3.1| - P7^ ). In virtue of ( p^.4lD and ( |3.15| ) 

T^A (3.17) 

where the constant A takes the value 10^ K^/^ m~^, Mo is the initial mass of star in meters, 
T the Kelvin temperature and p the dimensionless energy density. 



B. Maxwell-Cattaneo transport equations: explicit form 



The Maxwell-Cattaneo transport equations ( ^TIl - p.3|) can be written, by means of (1317 
and condition (|3.10| ), in dimensionless form 
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rf/°n „ + n ^ --aM^T^T 



'dP 



n 2 



15^ 



^lUl 



(3.18) 



(3.19) 



(3.20) 



where a ~ 6.252 x lO'^^ cm'^ K~^. 

The expression of mean free time for neutrinos ( p.l7|) is given for a comoving observer. 
Thus, applying a Lorentz boost to it, it is possible to write the transport equations for an 
observer with radial velocity —uo respect to matter, in Bondi coordinates. 

As we saw in last section, by means of the method HJR it is possible to solve the 
Einstein equations. Thus, the quantities Q, p, P^, -P± and u and the metric functions (3 and 
m are known. Resorting to expressions ( p.l7[ ), ( |2.1| ), ( |2.21| ) and ( |2.23| ) it is possible, after 
a straightforward calculation, to write down the Maxwell-Cattaneo equations in terms of 
known quantities Q, m, /?, tu, p and Pr and their derivatives. 
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The evolution of the heat flow is governed by 



f{u,r)i—^\ ~0o + /i(M,r)0i + M;(M,r) 



(3.21) 



where we have introduced for sake of simplicity 



and 
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(3.24) 
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(3.28) 



The bulk viscous pressure ( |3.2| ) can be found by solving the partial differential equation 
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(3.29) 
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where the expansion 9 can be written as 



m mi 



2m y/T^^uj^ \r'^ r 



2m 



+ e- 



-2/3 _ 



mo 



'1 -CJ 



l-LJ 



1 - 



2m 



(3.30) 



The description of shear viscous pressure evolution ( |3.3| ) is given, as in the other two 
dissipative flows, by a first order partial differential equation 
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where cr^ = |o"^i/Cr'^'^ and. 
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The transport coefficients (|3.7| - |3.9|) as measured by a comoving observer read. 



Pv04 



(3.31) 



(3.32) 



(3.33) 
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(3.35) 



with T ~ 4.173 x 10^^ g cm'^ K^/^ s. 

The system of three partial differential equations ( p. 211 ), ( p.29| ) and ( p.31| ) can be solved 
provided initials, finals and boundary conditions are given. The boundary condition for 
temperature in the heat conduction equation ( |3.21| ) deserves special attention. As shown in 
section IV C, the boundary condition for the temperature must fulfill some requirements. 
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After solving the heat conduction equation, we shall find the temperature, and using it in 



the evolution equations ( |3 .291 , 131311 ) the bulk and shear viscous pressure will be determined. 
Once the temperature is known, it will be easy to find the coefficients of heat conductivity, 
shear viscosity and bulk viscosity from equations (|3.33[-p.35D, respectively. 



IV. APPLICATION TO THE GOKHROO & MEHRA'S CONFIGURATION. 

As we have seen in section III, the H JR method starts from a static solution of Einstein's 
equations. It is possible to find numerical models based in nuclear physics as solutions of 
the Einstein equations. However, the uncertainties introduced by this procedure in the HJR 
method makes it desirable to start from a static analytical model instead. On the other 
hand, it is suitable to adopt an initially static fiuid distribution not excessively idealized. 
Some analytic anisotropic static solutions of the field equations are known |I|PJ39|-^. We 
shall adopt the static form of Gokhroo and Mehra (GM) |^ to illustrate the application of 
method described above. This solution corresponds to a fiuid with variable energy density 
in which some anisotropy is initially present even in absence of radiation. It leads, under 
some circumstances, to densities and pressures similar to the Bethe-Borner-Sato Newtonian 
state equation (BBS) H Jl0|j35|j4^ . Thus, this solution presents a compromise between the 



reahsm of the solutions based in nuclear physics and the desirable analytic expression of the 
static solution. 

A. Introduction of the GM configuration in the HJR method 

Following the first point of the HJR method (section II E), it is necessary to start from 
a static solution (p^t and P^t) of the Einstein equations. In the static solution of Gokhroo 
and Mehra p4| the energy density and radial pressure are assumed to be 

p(r) =pjl- k'pi (4.1) 



and 
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(4.2) 



where < k < 1 and n > 1 are constants. The central energy density and radial pressure 
are denoted by pc and Pc, respectively, and in the static case they are related by means of 
a constant A through the expression 



Thus, we identify pst as p{r) and Pgt as Pr{r). The tangential pressure is 

'p + Pr\ ( m{r) + Attv^PA 
^ r — 2m(r) J 

„2 



(4.3) 
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(4.5) 



with 



a — 



Snpc 



(4.6) 



In applying the HJR method (second point of the algorithm) we assume that the effective 
variables p{u,r) and P{u,r) have the same r-dependence than the physical quantities p{r) 
and P{r). The time dependence of these ones is introduced by two arbitrary functions of 
time K{u) and G{u). The form proposed here is 



p^ Pc{u) 1-K(u)^ , 



(4.7) 



and 



(4.8) 
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where 



Pc[U) = Pc 



Pc{u) = Pc 



K{u) 
K{u) 



K 



Pr 



K 



(4.9) 



(4.10) 



The central energy density and radial pressure in the final equilibrium state must differ from 
their values in the initial static case. Thus, we introduce their effective variables pc and Pc- 
These variables do not enjoy any physical meaning during the collapse, but they coincide 
with the central energy density and radial pressure in the static case. Because of this, it is 
necessary to introduce a time dependence in and P^.. The correctness of expressions ( [4.7| ) 



and (|4.8| ) may be seen a posteriori at the sight of the obtained results -section V below. As 



mentioned in subsection II E, the junction condition (|2.39|) allows us to relate both unknown 



functions of time K and G. Using (4/7) and ( [4.8|) in ( 2.39 ) and by means of ( 2.41 ) they are 
related by the expression 

(l-fi)(l-i^) 



(1 - GY 



FflX 



(4.11) 



To find an expression for m{u,r) and (3{u,r) (third point) we resort to the expressions 
( p.32| ) and (|2.35| ). In virtue of ( [4.7|) the mass function can be written as 
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where 



a 



K 



a- 



(4.13) 
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and the radius of the star a is written from ( |2.41| ) as A because, without loss of generality, 
the initial mass is taken as unity. The expression of P{u,r), by means of (|4.7| ) and ( [4. 81) , 
reads 
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(4.14) 
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where X may take different values depending on the relation between a. and K/A . 



12K ^ 
X takes the form 
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X takes a more simple form 

X=- 
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(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



Following the fourth point of the HJR method it is necessary to relate the unknown 
functions of time K and G with the functions A, Q, F or L. To this end, we evaluate the 
expression ([4.12|) in the surface of the star, 



K (aA^\ 



M = ma 

Ko \ 2 J \ 5 J 
On the other hand in virtue of ( p.41| ) the total mass of the system is 

M=|(l-F). 



(4.20) 



(4.21) 
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Equating the two last expressions it is possible to find the dependence oi K on A and F, 

(4.22) 



K{u)=K{A,F) = - 



1±1 1- 



fl2Ko\ fl-F 



\ 5a J \ 

The dependence of G on time (or equivalently on A, F and Q) can be established by 
substituting the last expression in ([4.11|) . Thus, introducing the expressions K{A,F) and 
G{A,F,n) in (|4^ ), (gj), (|4.12| ) and (|4.14| ), p, P, m and (3 are known throughout the star 
and for all times as functions of r and A, F and Q. 

To ascertain the physically meaningful sign in ( |4.22| ), we evaluate it at the initial time, 



(4.23) 



where we have made use of ( 2.41| ) with m(0) = 1. Therefore, if the initial value of the 



function K exceeds 5/6 we must adopt the positive sign in ([4.22|) , and the negative one 
otherwise. 

The functions A, F and Q can be found from the system of surface equations if the 



luminosity is known (point five of the method). The first two surface equations ( |2.42|) and 
( p.48| ), are model independent. In order to find a valid expression for the third one we resort 
to expression given in the last point of the algorithm (|2.53|) which in virtue of ([4.7| ), ( |4.8|) 



and ( |4.12| ) can be written as 
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Therefore, the system of surface equations is, 

A= Fin-1), 
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(4.25) 



(4.26) 



F= j[2L + F{1 - F){Q - 1)] 



(4.27) 
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and 



(4.28) 



where 
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(4.30) 



This system of differential equations can be solved for a given set of initial values of A, F 
and Q imposing a boundary condition (in this case the luminosity L). 

The initial central energy density can be written, by means of (|4.6| ) and ( [4.23|) , as a 



function of initial radius and Kn 



Pc 
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(4.31) 



AttAI (5 - 3Ko) 

Thus, a large radius implies a low central energy density. Moreover, the positive character 
of the discriminant in ( |4.22| ) implies some restrictions about the evolution of the total mass 
M of the star. During the collapse it must fulfill the inequality 

5a 



M < 



2ak; 



-A\ 



(4.32) 



In the static case a non negative value of K {i.e. Ko), ensures the condition dp/dr < 0, 
while condition Kq < 1 leads to a positive energy density. The effective energy density lacks 
of physical meaning when the star is collapsing. Therefore, the last argument cannot be 
applied to decide the range of values of function K here. Evaluating (|2.33|) in the surface 



and using ( p.40|) and ( |2.41| ) we obtain 



pa 



2Q- 1 



{Pa - Pro) ■ 



(4.33) 
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Q lies in the range 1/2 < Q < oo, further the condition p > P,. must be fulfilled. Thus, 
from the last expression, Pa > 0. Applying this restriction to ( |4.7| ) evaluated in the surface 
we obtain 



K{l-K)>{]^{)<K <l, 



(4.34) 



therefore the range of values physically admissible for K coincides with the static one. 
A last remark is in order before numerically solving the system of surface equations. Note 



that the generalization adopted for the tangential pressure ( p. 531) does not allow explosive 



models [34 . 



B. Evolution of the GM configuration 

We now are prepared to study the evolution of the initial static fiuid distribution de- 
scribed above in a particular case. To this end, it is necessary to adopt a specific initial 
configuration. The neutrinos can be studied in the diffusive regime if the density is higher 
than 10^^ ~ 10^^ g cm~^. As mentioned in the introduction, the neutrino trapping is an 
important source of viscosity. Because of that, in our diffusive model, the energy density in 
the surface must be, at least, about 10^^ g cm"'^. On the other hand a neutron star with a 
radius of about 10 km, has a central density not far from 10^^ g cm~^ [11|10|. In this diffusive 
model, the energy density in the center is a thousand times the corresponding to the surface. 
Thus, from (|4.1| ), Ko = 0.999. From expression (|4.31|) one can find the initial radius for a 



given initial stellar mass. Introducing the usual dimensions in (|4.31|) we obtain 

~ G47iAl{5-3K)Mr ^ ^^^^ 

where the initial mass Mo is given in geometrized units. Assuming Ao = 6, the initial mass 
is close to 1.3 M0. The corresponding initial radius can be found from ( p.41[ ) 



tta = AoMo ~ 11521 m. (4.36) 
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This values for the initial radius and mass are in the range of the usually accepted as typical 
for neutron stars. So we specifically adopt as initial configuration 



Pc ^ 1.01 X 10^^ g cm-^ 



Pa ~ 1.01 X 10^^ g cm 



Mo = 1.3 Mq 



and 



tto ^ 11521 m, 



(4.37) 



corresponding to a dense neutron star in which the diffusion approximation for neutrinos 
holds. 

To completely determine the initial conditions it is necessary to give a relationship be- 
tween the energy density and pressure in the center of the star ( [4.3| ). Also the parameter 
n occurring in the expression for pressure in the static case ( [4.2| ) must be imposed. If we 
assume the center of the star as a highly relativistic Fermi gas, then A = 1/3. In addition 
we take n = 1 for simplicity. 



Therefore, the initial conditions for the system of surface equations (|4.26| - ^l28|) are 



2 

3' 



1, Ko = 0.999, Pc = -pc and n = 1. 



(4.38) 



These ones correspond to a neutron star, initially at rest, with a redshift {za = l/y/F — 1) 
close to 0.22. 

The boundary condition, necessary to solve the system of surface equations, is supplied 
by the luminosity. We assume a Gaussian pulse for L centered in m = Upeak, and width A, 



L = - M-- 



Mr 



: exp 



u — 



A 



(4.39) 
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Note that Mr is the total energy radiated by the star in the collapse. The values adopted 
for Mr, A and Upeak are O.OIMq, 15 and 150 respectively. We have imposed the energy 



conditions for imperfect fluids |^ and the restriction — 1 < < 1 
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It is worth to emphasize the behavior of this model close to surface of the star. The 
surface collapses more slowly than the adjacent layers (figurel). It may be traced to the 
profile of the pressure gradient in this region. As seen in figure 2, its absolute value is larger 
in the surface than in the adjacent layers, which find a lesser resistance to collapse. Similar 
behavior is observed for lower values of Kg. Nevertheless, for low Ko the variation in the 
pressure gradient is smaller than for large Ko and, consequently, the difference between the 
velocity of surface and adjacent layers is not too high. 

The value of Ko largely determines the radius of the star once the collapse is finished. 
The higher Ko, the lower final radius. This behavior can be related to the fact that values 
of Ko far from unity follows from an equation of state close to incompressibility. Thus, in 
these cases, the velocity of the surface is not very high and possible variations in the radius 
are restricted. The highest variations in energy density and radial pressure take place in the 
most external layers, which comprise approximately a thousand meters below the surface 
(figures 3 and 4). 

The evolution of the energy density in the surface is depicted in figure 5. For Ko = 0.9 
{pa{u = 0) ~ 8.78 X 10^^ g cm-3) it varies about 6.5% , while for Ko = 0.999 the final 
energy density in the surface gets nine times its initial value. This behavior, with that of 
the velocity, hints that the compression of the star is higher in the shells near the surface. 
This is likely due to the radiation generated in the collapse. As it is showed in figure 6 the 
heat flow is more intense in the inner shells than in the more external ones. 

Initially the viscous pressure vanishes due to the absence of radiation. This condition 
is also fulfilled in the equilibrium situation once the collapse is over. Its contribution to 
the tangential pressure during the collapse must be established by means of the transport 
equations. 
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C. Boundary and initial conditions for the transport equations 



We have solved the transport equation for the heat flow ( 3.21 ) for different boundary 
conditions. The interior solution found for the temperature ceases to be sensitive to the 
boundary condition about 500 ~ 700 m below the surface. This hints that the energy of 
the emergent neutrinos is not correlated with its energy beyond this external region and, 
consequently, most of the neutrinos escaping the star have been generated in the inner limit 



of this region. This agrees with the idea of neutrinosphere P^. A possible cause of this 



behavior can be found in the abrupt decrease of the heat flow close to surface, that follows 
from the similar behavior shown by the energy density and radial pressure in this zone. As 
the model of Gokhroo and Mehra, the BBS model and others based on nuclear physics show 
the same features HJTO|]. Thus, the neutrinosphere seems to occur whenever the adopted 
model is not highly idealized. 

Aimed to establish the boundary condition for the temperature we introduce the effective 
temperature, Te//. This one is usually defined by means of 



E 



(4.40) 



where a = be /A, and 



E 



inr^e] _ (4.41) 



is the luminosity as measured by the non-comoving observer momentarily located on the 
surface. It is to say, Tg// would be the temperature in the surface of the star if it would 
radiate as a black body. The idea of effective temperature -see for instance page 586] 
and [^, page 295]- is applicable to the most external layers of the star. It is related with 



the material temperature by the expression 



= Ir^ff (l + ) ' (4-42) 



where r is the optical depth {i.e. dr = —dr/Xf,ff). 
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According to last expression, if r = 2/3 the effective temperature coincides with the 
material one. Thus, most of the emergent neutrinos are generated in a shell close to r ~ 2/3. 
This one is, in the model under consideration, about 500 meters below the surface. At the 
surface r vanishes. Hence, the material temperature of the surface can be written as 



rp4 



_ 1 r 
a ~ 2 



rri4 

^eff 



(4.43) 



The boundary condition can be found from last expression together with ( |2.43| ), ( p. 471) , 
(g) and (1^. It reads 

2fi - 1 



t: = Qa- 



2a 



(4.44) 



This condition on the surface affects the evolution of the temperature in the layers close to 
the surface but not the evolution of the inner temperature (figure 7). 

To find explicitly the bulk and shear viscous pressure it is necessary to impose initial 
and boundary conditions. Something about the initial condition for the viscous pressure 
has been noted at the end of the last subsection. We assume that viscous processes appear 
because of the interaction between matter and neutrinos. Therefore, initially 11 and tt must 
vanish because of the absence of radiation. 

At the center of the star the isotropy condition {P± = Pr) is fulfilled. Thus, in virtue of 
( p.l9| ) the relation 



lim a 7r 

r^o y 2 







(4.45) 



must be satisfied. And as a consequence of the regularity condition [mj^.^^ = 0, at the center 
of the star the generation of neutrinos is forbidden. Therefore we assume that both viscous 
pressures, vanish at r = 0. Thus, the boundary conditions imposed on ( p. 291) and ( p.31|) are 



[UU, = [nU, = Vu, 



(4.46) 



and from ([4.45|) 



Jr=0 



Vu. 



(4.47) 
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D. Temperature and viscous pressure: main results 



We have used the Maxwell-Cattaneo transport equations instead of Eckart's since the 
latter ones violate relativistic causality. The temperature in the inner layers found by means 
of ( 3.21| ) is the same as the one obtained using (|3.4|) (i.e., Fourier's law). Nevertheless, in 



the neutrinosphere the differences between both solutions cannot be neglected, especially in 
the first and late stages of the collapse (figure 8). The main difference between Maxwell- 
Cattaneo and Eckart heat transport equations is the presence of the term r Q!" h'^, which 
introduces in ( ^.21 ), the term 



fM\^j , (4.48) 

while leads to 

(4.49) 

In the neutrinosphere and in the case under consideration one has, f{u,r) ~ 0.1 g{u,r) and 
B ~ 10^^, i.e. both terms are of about the same order when temperature approaches 10^^ K. 
Thus, the effects introduced by the presence of a relaxation term in heat transport equation 
are important there. 

The high value reached by the temperature (figure 7) validates the adopted initial and 
final conditions: the approximation that in the static case the temperature for the neutron 
star vanishes. Applying this condition with ( [4.44| ) to expression ( p. 211) allows us to find the 



evolution of the temperature. Once the temperature is known, the transport coefficients can 
be calculated with the help of (|3.33| - |3.35D . 

The bulk viscous coefficient, (, vanishes at the center of the star because of the state 
equation ([4.38|) adopted there. In the shells close to the center the state equation is near 
to Pr = p/3. Thus, the value of ( is less than corresponding to shear viscous coefficient (rj) 
(figures 9 and 10). Nevertheless, in virtue of ( |3.9|) , in stars with a central equation of state 
different from the adopted in this case -i.e. not highly relativistic-, ( can be of the order of 
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1]. On the other hand a comparison among the terms in the evolution equations (|3.29|) and 
( p.31| ) points out that the term 

2 

(4.50) 



1 _ dP\ 
3" 



is responsible for the main differences between the bulk and shear viscous pressures. 

In the region in which the energy density is roughly 10^^ g cm^^ we find for the shear 
viscosity coefficient values about 10^*^ dyne cm~^ s~^, which is 10^ times higher than the 
corresponding to the interactions among electrons, protons and neutrons at the same density 
p6|| . This underlines the importance of the neutrino trapping as a source of viscosity in 
the stellar collapse. Similarly the thermal conductivity coefficient is greater for neutrino 
scattering than the mentioned . 



The evolution of vr is shown in figure 11 for different layers of the star. The shear viscosity 
value rises swiftly form zero, at the center of the star, and becomes an important source of 
anisotropy in the core of the star (figure 12). For distances from the center larger than 2 
km, vr can be neglected as compared with p as a source of anisotropy. 

The bulk viscous pressure is depicted in figure 13. In this model 11 can be neglected 
against vr in the innermost shells. However, in the peripheral layers 11 and vr, are of the same 
order. 

V. DISCUSSION 

In this paper we have used the HJR method to generate non static solution departing 



from the static anisotropic fluid distribution of Gokhroo and Mehra p4[. This model shows 
an initial non uniform energy density. The presence in the expression for the energy density 
of the parameter K allows us to model without difficulties the initial features of the star 
according to the radiation limit adopted on it (diffusion or free-streaming out). 

During the collapse huge quantities of neutrinos are generated. These ones transport 



thermal energy from the most interior regions to the exterior layers P,P|,pT]|. Because of the 
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high densities and temperature, the neutrinos interact with matter. In their trip toward 
the surface they get thermahzed and drive the star to a new equihbrium state. The most 
important thermahzation processes are the absorption of neutrinos by neutrons and colhsions 
between neutrinos and electrons. Consequently we have chosen a density range according to 
these interactions (pc ~ 10^^ g cm~^, pa ~ 10^^ g cm~'^), allowing us to treat the radiation 
in the diffusive limit. 

The density and pressure profile in the central layers vary little during the collapse. 
Therefore, the behavior of the zones nearby the surface is perhaps the most important. Due 
to the decrease, in these layers, of the pressure gradient with respect to the surface, the 
velocity of collapse of the former is larger than the corresponding to this one. This behavior 
also occurs in models with a density at the surface of about lO^^g cm~^, though then the effect 
is not so intense. Also, the final radius of the star depends on the density at the surface. The 
higher initial energy density at the surface, the larger final radius. As previously mentioned, 
the largest variations in energy density occur at the surface. This effect is probably due to 
the importance of the radiation in the evolution of the neutrinosphere. In that zone, the 
heat flow is much weaker than in the inner regions. Nevertheless, it is possible that the 
radiation energy density in it comprises a large part of total energy density. 

The knowledge of the temperature of the star is the key to establish which processes 
can take place in its bosom. At the variance with the heat flow, this quantity cannot 
be obtained from Einstein's equations. It is absolutely necessary then to appeal to some 
thermodynamic theory of irreversible processes. We have taken a further step with respect 



to other analysis by introducing the EIT in its covariant formulation [^,0, avoiding in 
this way the unphysical behavior of the conventional theory The rigorous study of 

transport equations derived from EIT is technically complex. Therefore, we have appealed 
to the so-called "truncated" form that reduces to the Maxwell-Cattaneo equations |2^. To 
solve the corresponding equation for temperature we have found the mean free time of the 
neutrinos. In so-doing we have considered both reactions mentioned above, if the neutrinos 
are generated by thermal emission with energy near ksT . In this way it is possible to 
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establish the dependence of the thermal conductivity coefficient on temperature (x oc T^/^). 
Then, solving the equation for the evolution of the heat flow, we found the temperature. 
Also we have solved the classic transport equation, finding out considerable differences in 
the neutrinosphere between both theories. These temperatures disagree from one another 
in the ffist and last stages of the collapse. 

It is worth to emphasize a peculiarity of the temperature in the innermost layers. This 
one is insensitive in that zone to the boundary condition imposed on the surface. Therefore, 
those neutrinos that manage to escape from the star lack information about the temperature 
of the central regions, in other words their energy in the surface is not correlated with that of 
the interior. This suggests that these have been created in a region close to the surface, the 
so-called neutrinosphere - [§, page 586], page 295], page 17]. The neutrinosphere is 
commonly associated with the effective temperature, Teff, which coincides with the material 
temperature in its lower limit. Through the introduction of T^ff we find out the temperature 
at the surface. Thus, we have obtained a credible boundary condition for the heat transport 
equation. We have supposed, therefore, that the layers of the star that are sensitive to a 
change in the boundary condition form the neutrinosphere with a thickness of 500 ~ 700 
meters. This length approximately coincides with the neutrino mean free path in that zone. 
This fact reinforces the hypothesis that the emergent neutrinos have been created in the 
interior region that delimits the neutrinosphere. 

Once the temperature was found, we went on by solving the transport equations of the 
remaining dissipative flows (bulk and shear viscous pressure). Using Weinberg expressions 
for the viscosity coefficients and expression ( p.l7|) , we have obtained C^V T^/^. This 
yields a value about 10^^ poise for shear viscosity coefficient. This value is much larger than 
the corresponding for interactions among electrons, protons and neutrons. Although ( can 
be of the same order that 1], the value we have obtained for ( is about one hundred times 
lower than rj. This is because we have assumed the center of the star as an ultrarelativistic 
Fermi gas and it implies an equation of state of the form P = p/3. It is well known that 
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for a fluid governed by this equation the bulk viscosity coeflicient vanishes 0. Though 
throughout the star this hmit is not strictly fulfllled, the relationship between the radial 
pressure and energy density restricts the value of In spite of this, from seven kilometers 
of the center upward the bulk viscous pressure is comparable to shear one. In an anisotropic 
model, shear viscous pressure is restricted by the difference between the tangential and radial 
pressures. Therefore, one can expect that its contribution to the total pressure is small. It is 
worth noting the contribution of shear viscous pressure to the anisotropy. According to the 
model studied here, it seems that the viscosity is responsible for an important part of the 
inner anisotropy. Then, the shear viscous pressure greatly contributes to total anisotropy 
in approximately the two nearest kilometers to the center of the star. It is certain that in 
the most internal zone of the star the anisotropy is less important that in the peripheral 
region. Though, under certain circumstances (larger collapse speeds) the importance of 
viscous pressure in connection to the anisotropy will be extended to more afar zones. 
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FIGURE CAPTIONS 

Figure 1. Velocity profile a for given constant Scliwarzscliild time (T = 0.9 ms) and two 
different values of Ko- All the figures made at constant Schwarzschild time were con- 
structed from expression (|2.7] ). 

Figure 2. Pressure gradient profile. The decrease of |Pr.i| close to surface is, probably, the 
cause of the behavior shown by the velocity in this region. Schwarzschild time T = 0.9 
ms. 

Figure 3. Relative variation of the energy density as a function of the retarded time u. 
The corresponding value for Ko is 0.999. po denotes the initial energy density. 

Figure 4. Relative variation of the radial pressure as a function of the retarded time u. 
Ko = 0.999. 

Figure 5. Energy density evolution at the surface for two different models {Ko = 0.999 and 
Ko = 0.9). Note the different scale used in each case. 

Figure 6. Radial heat fiow Q as measured by the Minkowskian observer locally comoving 
with matter. Ko = 0.999 and T = 0.95 ms. 

Figure 7. Temperature as measured by an observed momentanously located at the surface. 
Its maximum value at the surface is about 5.44 x 10^^ K. Here, and in the following 
figures, 3^e = 0.2. 

Figure 8. Difference in the neutrinosphere between the temperature found by means of the 
Eckart and Maxwell-Cattaneo transport equations. 

Figure 9. Evolution of the bulk viscosity coefficient (. 

Figure 10. Evolution of the shear viscous coefficient rj. 

Figure 11. Evolution of the shear viscous pressure in different layers of the star. 
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Figure 12. Comparative plot of tlie different contributions to the anisotropy (p.l9|) at r = 
O.lAo- Lines marked A and B correspond to terms —3n/2 and p, respectively. 

Figure 13. Evolution of the bulk viscous pressure within the star. 
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